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We present calculations of the quantum and thermal Casimir interaction between real mirrors in 
electromagnetic fields using the scattering approach. We begin with a pedagogical introduction of 
this approach in simple cases where the scattering is specular. We then discuss the more general case 
of stationary arbitrarily shaped mirrors and present in particular applications to two geometries of 
interest for experiments, that is corrugated plates and the plane-sphere geometry. The results nicely 
illustrate the rich correlations existing between material properties, temperature and geometry in 
the Casimir effect. 



I. INTRODUCTION 



versal expressions for the force -Fcas and energy Eq^ 
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The Casimir effect [l| is an observable effect of vac- 
uum fluctuations in the mesoscopic world, to be tested 
with the greatest care as a crucial prediction of quantum 
field theory d-Q. It also constitutes a fascinating inter- 
face between quantum field theory and other important 
aspects of fundamental physics, for example through its 
connection with the problem of vacuum energy [9-11]. 

Casimir physics plays an important role in the tests 
of gravity at sub- millimeter ranges 0, [l^l . Strong con- 
straints have been obtained in short range Cavendish-like 
experiments [ij] . A hypothetical new force of Yukawa- 
like form could not exceed the gravitational force in the 
range above 56//m. For ranges of the order of the mi- 
crometer, similar tests are performed by comparing the 
results of Casimir force measurements with theoretical 
predictions (T5l - [l7| . At even shorter scales, those tests 
can be performed using atomic [T^ or nuclear [T9| force 
measurements. In any of these short-range gravity tests, 
a new hypothetical force would appear as a difference 
between the experimental result -Foxp and the theoretical 
prediction i^th- This implies that i^th and -Foxp have to 
be assessed independently from each other and necessar- 
ily forbids to use the theory-experiment comparison for 
proving (or disproving) some specific experimental result 
or theoretical model. 

Finally, the Casimir force and the closely related Van 
der Waals force are dominant at micron or sub-micron 
distances, entailing their strong connections with vari- 
ous important domains, such as atomic and molecular 
physics, condensed matter and surface physics, chemical 
and biological physics, micro- and nano-technology pol |. 



II. COMPARISON OF CASIMIR FORCE 
MEASUREMENTS WITH THEORY 



Casimir calculated the force between a pair of perfectly 
smooth, flat and parallel plates in the limit of zero tem- 
perature and perfect reflection which led him to the uni- 
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with L the mirrors' separation, A their surface, c the 
speed of light and h the Planck constant. The univer- 
sality of these ideal Casimir formulas is explained by the 
saturation of the optical response of perfect mirrors which 
exactly reflect 100% of the incoming fields. This idealiza- 
tion does not correspond to any real mirror. In fact, the 
effect of imperfect reflection is large in most experiments, 
and a precise knowledge of its frequency dependence is 
essential for obtaining reliable theoretical predictions to 
be compared with Casimir force measurements j21i - t35i] . 



A. The description of metallic mirrors 

The most precise experiments are performed with 
metallic mirrors which are good reflectors at frequencies 
smaller than their plasma frequency wp. Their optical 
response at a frequency w is described by a reduced di- 
electric function written as 



(2) 



The function e [uS\ represents the contribution of inter- 
band transitions and it is regular at the limit w 0. 
Meanwhile a [uj] is the reduced conductivity, measured as 
a frequency (the SI conductivity is eoc), which describes 
the contribution of the conduction electrons. 

A simplifled description corresponds to the lossless 
limit 7 — often called the plasma model. As 7 is much 
smaller than wp for good conductors, this simple model 
captures the main effect of imperfect reflection. However 
it cannot be considered as an accurate description since 
a much better flt of tabulated optical data is obtained 
with a non null value of 7 (sgI . Wf\ . Furthermore, the 
Drude model, with 7 7^ 0, meets the important property 
of ordinary metals which have a flnite static conductivity 

(To = in contrast to the lossless limit which corre- 
sponds to an inflnite value for ctq. 
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Another correction to the Casimir expressions is as- 
sociated with the effect of thermal fluctuations (38l - l4ll |. 
Bostrom and Sernelius have remarked that the small non 
zero value of 7 had a significant effect on the force eval- 
uation at ambient temperature !42] . This significant dif- 
ference is attributed to the vanishing contribution of TE 
modes at zero frequency for dissipative mirrors entailing 
that for the Casimir force, contrary to the dielectric func- 
tion, there is no continuity from the Drude to the plasma 
model at the limit of a vanishing relaxation. The ratio 
between the predictions evaluated at 7 = and 7^0 
even reaches a factor 2 at the limit of large temperatures 
or large distances. Unfortunately it has not yet been 
possible to test this striking prediction since the current 
experiments do not explore this domain. 

The current status of Casimir experiments appears to 
favor theoretical predictions obtained with the lossless 
plasma model 7 = rather than those corresponding to 
the Drude model with 7 ^ as one might have expected 
(see Fig.l in [1^). We thus have to face a discrepancy 
between theory and experiment. This discrepancy may 
have various origins, in particular artefacts in the experi- 
ments or inaccuracies in the calculations. They may also 
come from yet unmastered differences between the situa- 
tions studied in theory and the experimental realizations. 

These remarks have led to a blossoming of papers de- 
voted to the thermal effect on the Casimir force, for re- 
views see e.g. (43l - l47j . It is worth emphasizing that mi- 
croscopic descriptions of the Casimir interaction between 
two metallic bulks lead to predictions agreeing with the 
lossy Drude model rather than the lossless plasma model 
at the limit of large temperatures or large distances [48l - 

It is also important to note that the Drude model leads 
to a negative contribution of the Casimir interaction to 
entropy, in contrast to the plasma model [5l| . There is no 
principle inconsistency with the laws of thermodynamics 
at this point since the negative contribution is^ nothing 
but a difference of entropies (see for example [52j). 



B. The role of geometry 



sphere geometry in particular, the PFA can only be valid 
when the radius R is much larger than the separation 
L [59I - I61I I. But even if this limit is met in experiments, 
the PFA gives no information about its accuracy for a 
given ratio of L/ R and how this accuracy depends on the 
properties of the mirror, on the distance or temperature. 

Answers to these questions can only be obtained by 
pushing the theory beyond the PFA, which has been done 
in the past few years |62h66| . A multipolar expansion of 
the Casimir effect between perfect mirrors in electromag- 
netic vacuum was proposed in [67l . [68l |. These calculations 
have now been performed for plane and spherical metallic 
surfaces coupled to electromagnetic vacuum, at zero 
or non zero temperature [t^. El , which has opened the 
way to a comparison with theory of the only experimen- 
tal study devoted to a test of PFA in the plane-sphere 
geometry [t^I ■ As we will see at the end of this article, 
the features of the thermal Casimir force mentioned in 
section |ll] are considerably altered when the geometry is 
properly taken into account. The factor of 2 between the 
force values within Drude and plasma model is reduced 
to a factor of 3/2 decreasing even more below this value 
when small spheres are considered. Negative entropies 
are not only found for the Drude model but also for per- 
fect reflector and plasma models, which means that neg- 
ative contributions of the Casimir interaction to entropy 
can be found even in the absence of dissipation. 

Another specific geometry of great interest, that we 
will present in the following, is that of surfaces with pe- 
riodic corrugations. As lateral translation symmetry is 
broken, the Casimir force contains a lateral component 
which is smaller than the normal one, but has neverthe- 
less been measured in dedicated experiments [73j . Calcu- 
lations beyond the PFA have first been performed with 
the simplifying assurnptions of perfect reflection [tJ] or 
shallow corrugations [tI, [t^ . As expected, the PFA was 
found to be accurate only in the limit of large corruga- 
tion wavelengths. Very recently, experiments have been 
able to probe the beyond-PFA regime [t^, Ii3 and ex- 
act calculations of the forces between real mirrors with 
deep corrugations [tI, H^I have been performed. More 
discussions on these topics will be presented below. 



The geometry plays an important role in the context of 
theory/experiment comparison for Casimir forces. Pre- 
cise experiments are indeed performed between a plane 
and a sphere whereas most exact calculations are devoted 
to the geometry of two parallel plates. The estimation 
of the force in the plane-sphere geometry thus involves 
the so-called Proximity Force Approximation (PFA) [53] 
which amounts to averaging the force calculated in the 
parallel-plates geometry over the distribution of local 
inter-plate distances, the force being deduced from the 
Lifshitz formula [sj, [5^, the meaning of which will be 
discussed below. 

This trivial treatment of geometry cannot reproduce 
the rich interconnection expected to take place between 
the Casimir effect and geometry [SGl - lssj . In the plane- 



Ill. THE SCATTERING APPROACH 

In the following, we will focus our attention onto the 
scattering approach, which is an efficient and elegant 
method for addressing the aforementioned questions. 

This method has been used for years for describing the 
optical properties of non perfectly reflecting mirrors in 
terms of scattering amplitudes [sil, [s^l- These scatter- 
ing amplitudes are often deduced from Fresnel reflection 
amplitudes calculated for mirrors described by local di- 
electric response functions, in which case the expression 
of the Casimir force is reduced to the Lifshitz expression 
(53 . [55I . However the scattering approach is much more 
general than the Lifshitz one since real mirrors are always 
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described by some scattering amplitudes but not neces- 
sarily by local dielectric response functions. This point 
will be discussed in more detail below. 

The interest in the scattering approach has consider- 
ably increased since it has become clear that it is also an 
extremely efficient method for calculating the Casimir 
effect in non trivial geometries. This was realized by 
several groups employing different theoretical techniques 
and using different notations (see for an historial 
overview). Besides the already quoted papers, one may 
cite the following references which used different versions 
of the scattering approach (s^ - fssj or alternative methods 
[sqI - Iq^ I . This topic has seen recently an impressive num- 
ber of new a pplications proposed, among ones one may 
cite [flMn^. 

The first explicit application of the scattering approach 
to non-trivia l geometr ies and non perfect refiectors was 
developed in [l03l Il04| to calculate the roughness correc- 
tion to the Casimir force between two planes, in a per- 
turbative expansion with respect to the roughness am- 
plitude. The same perturbative formalism was also ap- 
plied to compute t he la teral Casimir force [zl, 1?^ and 
the Casimir torque [fOSj between two corrugated surfaces 
made of real material, and then to derive the Casimir- 
Polder po tential for an atom near a corrugated surface 

[iMIinsj. 

Let us recall that results app li cabl e to the non retarded 
case have been available |lQ9l , Ill0j | before those corre- 
sponding to the full retarded theory, and also that the 
scattering theory has been used for a long time for study- 
ing the C asimir-Polder force between atoms or molecules 

[nillm 



We begin the review of the scattering approach by 
an introduction considering the two simple cases of the 
Casimir force between two scatterers on a 1-dimensional 
line and between two parallel plates coupled through 
specular scattering to 3-dimensional electromagnetic 
fields [81j . We then address the general case of non spec- 
ular scatterers in 3-dimensional electromagnetic fields Q . 



A. Mirrors on a 1-dimensional line 

The first case corresponds to the quantum field theory 
in 2-dimensional spacetime (1-d space plus time). In this 
simple case, we have to consider only two scalar fields 
counter-propagating along opposite directions. The re- 
sults summarized below are drawn from a series of papers 
devoted to the study of static or dynamic Casimir forc e 
betw een mirrors coupled to these scalar fiel ds M , fsil . Ill3l - 
Il2l| . For example, it was established in [ll8] that the 
Casimir energy does contribute to the inertia of the cav- 
ity as it should according to the principles of relativity. 

In this simple model, a mirror Mi at rest at position 
qi is described by a 2x2 scattering matrix Si containing 
reflection and transmission amplitudes ri and ti 

^ — 2iujqi / c 
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Two mirrors Mi and M2 at rest at positions qi and q2 
form a Fabry-Perot cavity described by a global scatter- 
ing matrix 6*12 which can be deduced from the elementary 
matrices Si and ^2 associated with the two mirrors. 
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The denominator d is given by 

d = 1 - rir2e2^"^/^ , L = q2 ^ qi, 



(5) 



and its zeros (the poles of S'12) represent the resonances 
of the cavity. It turns out that the forthcoming discus- 
sions of the Casimir effect depend only on the expression 
of d and not on all the other details in the form of S'12. 
The reason explaining this property is the following re- 
lation between the determinants of the S— matrices (all 
supposed to be unitary in the simple model) 



detSi2 = (detSi)(detS2) 



(6) 



From this relation, it is easy to derive the Casimir free 
energy as a variation of field energy (vacuum energy at 
T = 0, vacuum plus thermal energy otherwise). The 
presence of a scatterer indeed shifts the field modes and 
thus induces a variation of the global field energy. The 
Casimir free energy is then obtained as the variation of 
field energy in presence of the cavity corrected by the 
effects of each mirror taken separately 
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A is a function of the frequency to representing the phase- 
shift produced by the Fabry-Perot cavity, again corrected 
by the effects of each mirror taken separately 



A(^) 



lndetSi2 - IndetSi - lndetS2 



(8) 



N is the mean number of thermal photons per mode, 
given by the Planck law, augmented by the term ^ which 
represents the contribution of vacuum 
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This phaseshift formula can be given alternative inter- 
pretations 81]. In particular, when the Casimir force F 
is derived from the free energy 
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it is seen as resulting from the difference of radiation 
pressures exerted onto the inner and outer sides of the 
mirrors by the field fluctuations. For each field mode at 
frequency w, represents the field momentum while g 
is the ratio of fluctuation energies inside and outside the 
Fabry-Perot cavity. 

Using the analytic properties of the causal function 
In d, the Casimir free energy can also be written as an in- 
tegral over imaginary frequencies uj — (Wick rotation) 



J- = h I — cot ( —T-T^ 

2tt \2k^T 



(11) 



Using the pole decomposition of the cotangent function 
and the analytic properties of Ind, this expression can 
finally be written as a sum over Matsubara frequencies 
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The Matsubara sum is the sum over positive integers 
m with m — Q counted with a weight ^. 

For completeness, let us recall also that the contri- 
bution to entropy of the Casimir interaction is simply 
written as 



S = - 



dJ-{L,T) 
df ■ 



(13) 



Hence, it is defined as a difference of entropies just as the 
free energy T has been defined in ([7]) above as a difference 
of free energies. 



B. Specular reflection in 3-d space 

The same lines of reasoning can be followed when 
studying the case of two specularly reflecting mirrors 
coupled to electromagnetic fields in 3-dimensional space. 
The geometry is sketched in Fig. [T]with two plane paral- 
lel mirrors aligned along the transverse directions x and 
y (longitudinal direction denoted by z) 



FIG. 1: Two plane parallel plates at distance L facing each 
other constitute the Casimir cavity. 



Due to the symmetry of this configuration, the fre- 
quency w, the transverse vector k = {kx,ky) and the 
polarization p = TE, TM are preserved by all scattering 
processes. The mirrors are described by reflection and 
transmission amplitudes which depend on these scatter- 
ing parameters. We assume thermal equilibrium for the 
whole "cavity -I- fields" system, and proceed with the 
derivation as in the simpler case of a 1-dimensional space. 
Some elements have to be treated with greater care now 
[1, [s^l . First there is a contribution of evanescent waves 
besides that of ordinary modes freely propagating out- 
side and inside the cavity and it has to be taken carefully 
into account. The properties of the evanescent waves are 
described through an analytical continuation of those of 
ordinary ones, using the well defined analytic behavior 
of the scattering amplitudes. Then dissipation inside the 
mirrors may also play a role which implies considering 
the additional fluctuation lines coming along with dissi- 
pation ^8, 82]. 

At the end of this derivation the free energy may still 
be written as a Matsubara sum 



= fceT ^ ^ ^ ' In d{i^,n, K p) 



(14) 
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is the sum over transverse wavevectors with A the 
area of the plates, the sum over polarizations and 

' the same Matsubara sum as in the 1-d case. The 
denominator is now written in terms of the result k of 
Wick rotation on the longitudinal wavevector k^ 



dii^Kp) = 1 - 7-i(«C,k,_p)r2(iC,k,p)exp 



-2kL 



(15) 



k2 + ^. 



This expression reproduces the ideal Casimir formula 
([T]) in the limits of perfect reflection rir2 — J> 1 and zero 
temperature T — > 0. It is valid and regular at ther- 
mal equilibrium at any temperature and for any optical 
model of mirrors obeying causality and high frequency 
transparency properties d, [Sll, [13] . It can thus be used 
for calculating the Casimir force between arbitrary mir- 
rors, as soon as the reflection amplitudes are specified. 
These amplitudes are commonly deduced from models 
of mirrors, the simplest of which is the well known Lif- 
shitz model [54, 55] which corresponds to semi-infinite 
bulk mirrors characterized by a local dielectric response 
function £{uj) and reflection amplitudes deduced from the 
Fresnel law 



rTE(k,0 = 



K + Kt 



rTM (k,0 



en — Kt 



£K - 



,(16) 



(17) 



e is the dielectric function ^ and Kt denotes the result of 
Wick rotation of the longitudinal wavevector inside the 
medium. 
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In the most general case, the optical response of the 
mirrors cannot be described by a local dielectric response 
function. The expression (fT4)) of the free energy is still 
valid in this case with the reflection amplitudes to be 
determined from microscopic models of mirrors. Recent 
attempts in this direction can be found for example in 

At this stage, several remarks can be addressed to the 
readers interested in historical details: 

• The Lifshitz expression was not written in terms of 
reflection amplitudes un til K ats noticed that this 
formulation was natural (l29j |. To our best knowl- 
edge, the first appearance of an expression of the 
Casimir effect in terms of reflection amplitudes cor- 
responding to an arbitrary microscopic model (no t 
necessarily a dielectric response function) is in [8l| . 

• The fact that the expression of the free energy 
is valid for lossy as well as lossless mirrors is far 
from obvious. In the lossy case, one has indeed to 
take into account the contributions of fluctuations 
coming from the additional modes associated with 
dissipation. This property has been demonstrated 
with an incre asing range of validity in (8ll | , [s^ and 
[1] (see also [iSd ] for a theorem playing a crucial 
role in this demonstration). 

• The question had been asked in whether the 
regularity conditions needed to write the Matsub- 
ara sum were met for the Drude model which shows 
discontinuities at ^ — J' 0. This question has been 
answered positively in (52j |. 



tering matrix 



/>00 1 

= ih — iV(cj)lndetS 
Jo 27r 

ih I — iV(w)TrlnS. 
27r 



(18) 



The symbols det and Tr refer to determinant and trace 
over the modes of the scattering matrix at a given fre- 
quency u!. After a Wick rotation the formula can still be 
written as a Matsubara sum 

T=kBTJ2'TrlnD{^U , ^ (19) 



The matrix D (here written at Matsubara frequencies 
'-^m — i^m) is the denominator of the scattering matrix. 
It describes the resonance properties of the cavity formed 
by the two objects 1 and 2 and may be written as 



D 



1 - i?i exp"^^ i?2 exp"^^ . 



(20) 



The matrices i?i and R2 represent reflection on the two 
objects 1 and 2 respectively while exp~^^ describes 
propagation inbetween reflections on the two objects. 
Note that the matrices D, Ri and R2, which were diago- 
nal in the plane wave basis for specular scattering, are no 
longer diagonal in the general case of non specular scat- 
tering. The propagation factors remain diagonal in this 
basis with their eigenvalues k written as in (fH)) . Clearly 
the expression does not depend on the choice of a 
specific basis. We remark also that (fTO|) takes a simpler 
form in the limit of zero temperature 



dE 
dL 



E 



^Indei D{i(,). (21) 
27r 



C. The non-specular scattering formula 

We now present a more general scattering formula al- 
lowing one to calculate the Casimir force between sta- 
tionary objects with arbitrary shapes. We restrict our 
attention to the case of disjoint objects, exterior to each 
other, which corresponds to the configuration initially 
considered by Casimir (for interior configurations, which 
may be t reated with similar techniques, see for example 

HMH). 



The main generalization with respect to the already 
discussed specular cases is that the scattering matrix 
S has now to account for non-specular reflection. It 
is therefore a much larger matrix which mixes different 
wavevectors and polarizations while preserving frequency 
as long as the scatterers are stationary Q . Of course, the 
non-specular scattering formula is the generic one while 
specular reflection can only be an idealization. 

As previously, the Casimir free energy can be written 
as the sum of all the phaseshifts contained in the scat- 



Applications to be presented in the next sections will 
also involve the Casimir force gradient G which is often 
measured in experiments and defined as 



G=- 



dF 



(22) 



A number of the following applications will be discussed 
within the zero temperature limit, with a change of no- 
tation from the free energy !F to the ordinary energy E 
at zero temperature. 



IV. APPLICATIONS TO NON TRIVIAL 
GEOMETRIES 

Formula (PT|) has been used to evaluate the effect of 
roughness or corrugation of the s urfaces on the value of 
the Casimir force [Sll, [zl, [tI, Il03l | in a perturbative man- 
ner with respect to the roughness or corrugation ampli- 
tudes. It has also allowed one to study a Bose-Einstein 
condensate use d as a local probe of vacuum above a nano- 
grooved plate (lOTl Il08( . The scattering approach has 
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clearly a larger domain of applicability, not limited to the 
perturbative regime, as soon as techniques are available 
for c omp uti ng th e large matrices involved in its evalua- 
tion [ralsolfissj. 

Another important application, which we will summa- 
rize also in the present section, corresponds to the plane- 
sphere geometry used in most Casimir force experiments 
and for which explicit "exact calculations" (see a dis- 
cussion of the meaning of this expression below) have 
recently become available j67l - [7ll |. 



A. Perturbative treatment of shallow corrugations 

As already stated, the lateral Casimir force between 
corrugated plates is a topic of particular interest. It could 
in particular allow for a new test of Quantum ElectroDy- 
namics through the dependence of the lateral force on 
the corrugation wavevector (tsI . [76j . 




The part of the Casimir energy able to produce a lat- 
eral force is then found to be 



^lat 
J ^corrug 



dSE" 



db 



(24) 



KL 



■SR. 



exp" 



-KL 
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5Ri and 5R2 are the first-order variation of the reflection 
matrices Ri and R2 induced by the corrugations; -Dpiano 
is the matrix D evaluated at zeroth order in the corru- 
gations; it is diagonal on the basis of plane waves and 
commutes with K. 

Explicit calculations of ([M]) have been performed for 
the simplest case of experimental interest, with two cor- 
rugated metallic plates described by the plasma dielectric 
function. These calculations have led to the following ex- 
pression of the lateral part of the Casimir energy 

^^conug ^ ^Gc{kc)aia2 cos{kcb). (25) 

The spectral sensitivity function Gc(fec) has been given 
and discussed in (t!]. Using its expression, it is possible 
to prove a properly defined "Proximity Force Theorem" 
which states that the PFA is recovered at the limit of 
long corrugation wavelengths fee ^ 0. Obviously, this 
theorem does not imply that the PFA is always valid or, 
in other words that Gc(fcc) may be replaced by Gc(0). 



FIG. 2: Parallel corrugated surfaces, with L representing the 
mean separation distance, ai and 02 the corrugation ampli- 
tudes and b the lateral mismatch between the crests. When 
the corrugation are supposed to be the smallest length scales, 
the effect of the corrugations can be studied in the pertur- 
bative expansion. This approximation will be dropped later 



Here, we consider a geometry with two plane mirrors. 
Ml and M2, having corrugated surfaces described by uni- 
axial sinusoidal profiles such as shown in Fig. [21 We 
denote hi and /12 the local heights with respect to mean 
planes zi = and Z2 — L 



hi — ai cos(fccx) , /i2 
fee = 27r/Ac . 



a2 cos {kc{x - b)) , (23) 



hi and /12 have null spatial averages and L is the mean 
distance between the two surfaces; hi and /i2 are both 
counted as positive when they correspond to a decrease in 
the separation; Ac is the corrugation wavelength, kc the 
corresponding wavevector, and b the spatial mismatch 
between the corrugation crests. At lowest order in the 
corrugation amplitudes, when 01,02 <C Ac,Ap,L (with 
Ap the plasma wavelength describing the properties of the 
metallic mirror) , the Casimir energy may be obtained by 
expanding up to second order the general formula (EIJ. 
This perturbative approximation will be dropped in the 
next subsection. 



CI. 0.5 




FIG. 3: Variation of p versus the dimensionless variable 
kcL for metallic mirrors described by the plasma model, for 
fcpL =1 (dashed line), 2.5 (dotted line), 5 (dashed-dotted line) 
and 10 (solid line) [colors online with respectively green, blue, 
red and black lines]. 

To assess the validity of the PFA for the lateral Casimir 
force description, we now introduce the dimensionless 
quantity 



pike) 



Gc(fcc) 
Gc(0) ■ 



(26) 
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The variation of this ratio p with the parameters of in- 
terest is shown in Fig. [3] for gold covered plates with 
Ap = 137nm. The ratio p is smaller than unity as soon 
as kc significantly deviates from which means that the 
PFA overestimates the lateral Casimir force. For large 
values of fccj it even decays exponentially to zero, lead- 
ing to an extreme deviation of the real lateral force from 
the PFA prediction. 

Another situation of interest arises when the corruga- 
tion plates are rotated with respect to each other. As- 
suming as previously corrugations of sinusoidal shape 
with corrugation wavevectors k-,- having the same mod- 
ulus k — 27r/A(7 on both plates, it is possible to derive 
the second-order correction §E^™^^^ to the Casimir en- 
ergy which depends on the angle 9 between the corruga- 
tion s and thus has the ability to induce a Casimir torque 
[tsI llOSl j. Only crossed terms, proportional to the cor- 
rugation amplitudes on both plates, contribute to this 
expression, as the square terms are independent of the 
angle 9. The expression 5E'^°'^^^'' contains as the special 
case 9 — Q the pure lateral energy discussed above. Note 
that the dependence on the material properties and cor- 
rugation wavevector are captured by the same response 
function Gq already calculated. 

For quantitative estimations, we assume that the cor- 
rugations are restricted to a rectangular section of area 
LxLy with transverse dimensions and Ly much larger 
than the plate separation L and neglect diffraction at the 
borders of the plates. In Fig. SI we plot (5iJ'°'''5"° obtained 
in this manner, in arbitrary units, as a function of b and 
9. The Casimir energy is found to be minimal at 6 = 
and & = 0, Ac, 2Ac, which corresponds to a situation 
where corrugations are aligned. Starting from 9 — b — 
and rotating plate 2 around its center, one follows the 
line 6 = in Fig. H) Clearly, for small angles the plate is 
attracted back to 9 = b — without sliding laterally. 




-Ac 

FIG. 4: Casimir energy (arbitrary unit) as a function of the 
rotation angle 6 and the lateral displacement b. 

The Casimir torque is then deduced by deriving the 
energy with respect to the angle 9 

r = ^^6E'°"i'"'. (27) 



It is maximum at 9 = 0.66Xc / Ly where it is given by 

= 0.109 aia2kGcik)Ly. (28) 

The maximum torque per unit area is proportional to the 
length Ly of the corrugation lines, which plays the role 
of the moment arm. 

In contrast with the similar t orqu e appearing between 
misaligned birefringent plates [l36|, the torque is here 
coupled to the lateral force. This could induce com- 
plicated behaviours in an experiment and would prob- 
ably have to be controlled. This can be clearly seen 
on Fig. 01 if the plate is released after a rotation of 
9 > Ac /Ly it will move in a combination of rotation and 
lateral displacement. The energy correction vanishes at 
^ = ^c/Ly, defining the range of stability of the con- 
figuration b = 9 = 0. Rotation is favored over lateral 
displacements only for 9 < Xc/Ly. 




100 1000 
L (nm) 



FIG. 5: Maximum torque per unit area as a function of 
the mean separation L for the following parameters: aia2 = 
200 nm^, Ly = 24 ^im, Ap = 137 nm. Solid line: Ac = 2.4 ^m; 
dashed line: Ac = 1.2 /im; dotted line: Ac = 2nL/2.G (corre- 
sponding to the optimum value). 

However, the advantage of the configuration with cor- 
rugated plates is that the torque has a larger magnitude. 
Fig. [S] shows the maximum torque as a function of mean 
separation between the two corrugated gold plates with 
a plasma wavelength Ap = 137nni. At a plate separa- 
tion of about lOOnm the torque per unit area can be as 
high as 10~''Nm~^ These results on lateral forces and 
Casimir torques suggest that non trivial effects of geom- 
etry, i.e. effects beyond the PFA, can be observed with 
dedicated experiments. It is however difficult to achieve 
this goal with corrugation amplitudes oi, 02 meeting the 
conditions of validity of the perturbative expansion. This 
approximation is dropped in the next subsection. 
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B. Non perturbative calculations with deep 
gratings 

As already stated, recent experiments have been able 
to probe the beyond-PFA regime with deep corrugations 
[TTI Its! and it has also become possible to calculate exact 
expressions of the forces between nanostructures with- 
out using the perturbative assumption. This necessar- 
ily involves the non specular scattering formula and 
the evaluation of scattering properties mixing different 
wavevectors and polarizations. 

In the following we briefly discuss the Casimir interac- 
tion energy in a typical device made of two nanostruc- 
tured surfaces of intrinsic Silicon, such as shown in Fig. 
1 



did 



FIG. 6: Two surfaces with rectangluar gratings of depth h, 
gap width d and trench width d — di. 

To model the material properties of intrinsic Silicon, 
we use a Drude-Lorentz model f or w hich the dielectric 
function is well approximated by [l37l [ 



(29) 



with Eq k, 11.87 the value of the dielectric function at zero 
frequency, Eoo ~ 1.035 the high frequency limit of the di- 
electric function and wq = ^Co ~ 4.34 eV. Calculated with 
the proximity force approximation, the Casimir force be- 
tween the two gratings is given by the geometric sum of 
two contributions corresponding to the Casimir force be- 
tween two plates Fpp at distances L and L — 2h, which 
is independent of the corrugation period d. 

To assess quantitatively the validity of PFA, we plot 
as before the dimensionless quantity 



F 



FpFA 



(30) 



Fig. [7] shows this ratio for two Silicon gratings, separated 
by L = 250nm, of height h = lOOnm as a function of the 
corrugation period d with di = d/2 [7§|. Clearly, PFA is 
not a valid approximation except for two limiting cases, 
that is a vanishing corrugation period d — > and a very 
large corrugation periods d — >■ cxd, meaning in either case 
that the structured surfaces becomes flat. In between the 




20 40 60 



100 120 140 



FIG. 7: Casimir force normalized by its PFA value for two 
gratings of intrinsic Silicon with amplitude h = lOOnm and 
di = I as a function of d at a fixed distance L = 250nm. 



exact result for the Casimir force is always smaller than 
the PFA prediction, meaning that PFA overestimates the 
force. This has to be contrasted with calculations for 
perfect conductors where PFA always underestimates the 
real force. 

One important parameter to keep in mind is the num- 
ber of diffraction orders that has to be retained in the 
calculation in order for the Casimir energy to converge in 
the numerical calculation. This is illustrated in Fig. [5] for 
two Silicon gratings. For the sake of convenience, we plot 
the Casimir energy normalized by the energy for perfectly 
reflecting plane mirrors i.e. the energy reduction factor. 
The blue curve corresponds to the situation of two grat- 
ings of period 400nm separated by a distance L = 50nm. 
Clearly around five orders of diffraction are sufficient for 
the calculation of the Casimir energy in this case. The 
number of necessary diffraction orders decreases with in- 
creasing distance between the gratings. This is illustrated 
by the red curve where the two aforementioned gratings 
are now separated by a distance L — 400nm and where 
the Casimir energy has basically converged to its final 
value with only one order of diffraction retained. The 
fast convergence is here due to that fact that oblique 
diffraction orders are exponentially suppressed with in- 
creasing distance [l08j |. Finally, the greater the period 
of the grating the more orders of diffractions are needed 
as shows the green curve where the period of the two 
gratings is now 2/im. In this case, the Casimir energy 
has not yet fully converged to its final value even with 
as much as 13 orders of diffraction. This can be under- 
stood because the momentum transferred by the grating 
q = ^ is now small so that different orders of diffraction 
are nearly collinear with the specular one and therefore 
greatly contribute to the final energy. 

If attention is paid to the issue of convergence this cal- 
culation method is essentially exact and allows for direct 
comparisons with experimental results. In a recent ex- 
periment, Chan et al. have measured the Casimir force 
gradient between a gold sphere and a grating of doped sil- 
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FIG. 8: Convergence of the calculated Casimir energy be- 
tween two gratings as a function of the number of diffraction 
orders retained in the calculation. Gratings with different pe- 
riods are plotted as blue and red (400nm) and green (2 /im) 
points. The convergence of the calculations becomes slower 
when increasing the grating period d or decreasing the sepa- 
ration L. 



FIG. 9: Comparisons between experimental measurements 
and exact calculations for the Casimir force gradient between 
a gold sphere and two types of silicon gratings. Green and red 
dots correspond to data points provided by Ho Bun Chan for 
a grating period of 400nm and 1/im respectively. The solid 
curves of the same color are calculated data obtained using 
the scattering approach for the corresponding experimental 
parameters. 



icon [77|. Two samples of silicon gratings have been used. 
Both have a corrugation depth of l/xm, but different pe- 
riods of 400nni and 1/im respectively. The experimental 
data points of the ratio between the force gradient and 
its PFA approximation for both samples have been kindly 
provided by Ho Bun Chan and are plotted in Fig. [9l 

Concerning the calculation we model the optical prop- 
erties of Silicon by the dielectric function (|29)) . We have 
also taken into account the doping of the Silicon by 
adding a Drude part to this dielectric function, but this 
has not led to noticeable changes for the Casimir inter- 
action in the distance range up to 500nm which has been 
explored in the experiment. To model the optical proper- 
ties of gold we have used available optical data, extrapo- 

lated at low frequencies by a Drude model e(i^) = 
with ujp = 9eV and 7 = 35meV. The method is described 
in detail in [s^. The calculations were run up to = 3 
diffraction orders, after which the result for the Casimir 
energy was found to have converged. The result of our 
calculation is given in Fig. [S] as the solid green and red 
curves for the 400nm and l^m samples respectively. The 
theoretical predictions and the experimental data points 
are in good agreement. Due to a new improved numeri- 
cal code the agreement is better than the one presented 
in M. 



ing mirrors [671 168| where it was found that the Casimir 
energy was smaller than expected from the PFA and, 
furthermore, that the result for electromagnetic fields 
was departing from PFA more rapidly than was expected 
from previously existing scalar calculations [63 . [65l | . It is 
only very recently that the same calculations have been 
performed for the more realistic case of metallic mirrors 
at zero temperature [69l | and at arbitrary temperature 
[70I [71I I where both the lossless plasma model dielec- 
tric function and the lossy Drude dielectric function have 
been studied. We will sketch the method in the following. 




C. Exact calculations in the plane-sphere geometry 



FIG. 10: The geometry of a sphere of radius R and a plate at 
distance L; the center-to-plate distance is £, = L + R. 



The plane-sphere geometry is the configuration in 
which the most precise Casimir force measurements are 
currently performed [t^ . The Casimir interaction in this 
geometry can also be calculated in a formally exact man- 
ner using the general scattering formula (flQl) . Such cal- 
culations have first been performed for perfectly refiect- 



The set-up of a sphere of radius R above a flat plate 
is schematically presented in Fig. [TUl We denote respec- 
tively L and C = L + R the closest approach distance 
and the center-to-plate distance. In this configuration, 
the general expression of the Casimir free energy at tem- 
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perature T may be written as 

ra 

D = l- Rse-^^Rpe~^^. 



(31) 



The expression contains the reflection operators of the 
sphere i?s and the plate Rp which are evaluated with 
reference points placed at the sphere center and at its 
projection on the plane of the plate. They are sandwiched 
in between operators e~^^ describing the propagation 
between the two reference points. 

The upper expression is conveniently written through 
a dec omp osition on suitable plane-wave and multipole 
basis w]; Rp is thus expressed in terms of the Fresnel 
reflection coefficients rp with p =TE and TM for the two 
electromagnetic polarizations, while i?s contains the Mie 
coefficients , bi for respectively electric and magnetic 
multipoles at order £ — 1,2, .... Due to rotational symme- 
try around the z-axis, each eigenvalue of the angular mo- 
mentum m gives a separate contribution to the Casimir 
free energy obtained through the same formula as 

(I3ip . The scattering formula is obtained by writing also 
transformation formulas from the plane waves basis to 
the spherical waves basis and conversely. 

The result takes the form of a multipolar expansion 
with spherical waves labeled by £ and m {\m\ < £). It 
can be considered as an "exact" multipolar series of the 
Casimir free energy. Of course, the numerical computa- 
tions of this series can only be done after truncating the 
vector space at some maximum value f max of the orbital 
index i. 

The effect of this truncation is represented on Fig. [TT] 
where the Casimir energy in the plane-sphere geometry 
divided by its PFA estimation 



E 



PE 



^PFA 



(32) 



is plotted for various values of ^max, vly the special case 
of perfect mirrors in vacuum (T = 0). The figure shows 
that as expected the numerical results are more and more 
accurate when fmax is increased. More precisely the ac- 
curacy is significantly degraded when the ratio L/ R goes 
below a minimal value inversely proportional to ^max 



R 



1 



(33) 



To illustrate the effect of the truncation, one may say 
that the accuracy is degraded by typically more than 
0.1% when x < 0.05 for a value of ^max =85. For small 
values of x, which corresponds to the most precise current 
experiments, it may be possible to obtain information 
through an extrapolation of the numerical results. As an 
example, the dashed line on Fig. [TT] shows the result of 
a third degree polynomial fit using accurate numerical 
evaluations. 

As a further step, we show now on Fig. [12] the results 
corresponding to perfect and plasma mirrors, still at zero 
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FIG. 11: Upper graph: the ratio pE ~ E/E^^^ of the plane- 
sphere Casimir energy to its PFA estimation is plotted as a 
function of ^max for different values of L/R — 0.05,0.1,0.2. 
Lower graph: same ratio pE plotted as function of L/R for 
different values of fmax = 20, 40, 80. 



temperature [69|. We have derived the Casimir energy 
(jSTl) to obtain expressions for the force F and force gra- 
dient G, and then formed the ratios of the plane-sphere 
exact results to the PFA expectations F^^^ and G^^^ 
respectively 



PF 



pPFA 



PG 



G 



QPFA ■ 



(34) 
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FIG. 12: Variation of pa as a function of L/R as a function of 
L/R, for a nanosphere of radius R = lOOnm; the solid green 
line corresponds to gold-covered plates (Ap — 136nm) and 
the dashed red line to perfect reflectors. The decrease at low 
values of L/R represent a numerical inaccuracy due to the 
limited value of ^max (24 in this calculation 69|). 

Using these theoretical evaluations, it is now possible 
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to extract information of interest for a comparison with 
the experimental study of PFA in the plane-sphere geom- 
etry [22] ■ In this experiment, the force gradient has been 
measured for various radii of the sphere and no deviation 
of PFA was observed. The authors expressed their result 
as a constraint on the slope at origin /3g of the function 
Pg{x) 



Pg{x) = 1 + I3gx + 0{x'^) 



|/3g| <0.4. 



(35) 



Reasoning along the same lines, we have interpolated our 
theoretical evaluation of po at low values oix — L/R J3£] . 
Surprisingly the slope obtained for perfect reflectors was 
found to lie outside the experimental bound of [tI] 



-0.48. 



(36) 



The consistency with this bound is however recovered for 
the calculations done for plasma mirrors 



^plas 

Pn 



-0.21. 



(37) 



As a last example of application, we now discuss the 
effect of a non zero temperature. To this aim we evaluate 
eqn pri) at ambient temperature (T = 300K). The results 
of the numerical computations are shown on Fig. I13[ for 
the limiting case of perfect reflection (left) and for Drude 
metals (right) evaluated for Ap = 136nm, A-y/Ap = 250 
(values corresponding to gold). We have calculated the 
Casimir force F'^°^^ and F^rud j-jg^-y^reen the plane and 
the sphere at ambient temperature and then plotted the 
corresponding ratios 'ff?^'^^ and -jjOrud ^^^^^ force to a 
reference force corresponding to zero temperature 



Qjrpcrf 



erf 



pDrud 



(L,r) ^ 



dL 

gj-Drud 

dL 



Drud _ 



FP"^L,T) 

i^P«'-f(L,0) 

FD™d(x,r) 



(38) 



The various solid curves are drawn for different sphere 
radii i? as a function of the separation L ; the dashed 
curves on Fig. [T3] represent the quantities t^ppX ^^'^ ^pfa 
obtained from ((38)) by using the PFA ; the dotted curve 
in the upper graph is an analytical asymptotic expression 
discussed below. We do not show the corresponding plots 
for plasma mirrors as they are very similar to the perfect 
mirror case. 

The comparison of -(?p°''^ and ^^Diud j-gyeals surprising 
features, which could not be expected from an analysis in 
the parallel-plate geometry. First both ratios -d start from 
unity at small distances L/R 0. For R small enough, 
the ratios then decrease below unity with increasing dis- 
tance, reach a radius-dependent minimum and then in- 
crease again at very large distances. This behavior entails 
that the absolute value of the Casimir force is smaller at 
T = 300K than at T = 0, implying a repulsive contribu- 
tion of thermal fluctuations. The dashed PFA curve in 
the upper graph of Fig. [T3] representing t?pp^ is always 
larger than unity, excluding such a repulsive contribution 
of thermal fluctuations in the plane-plane geometry. 
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FIG. 13: Thermal Casimir force at T = 300K divided by the 
zero temperature force, computed between perfectly reflecting 
sphere and plane (upper graph), and between Drude metals 
(lower graph) plotted for Ap = 136nm, A^/Ap = 250. The 
solid lines from bottom to top correspond to increasing values 
of sphere radii. The dotted curve in the upper graph is the 
analytical asymptotic expression in the L S> -R limit. The 
PFA expressions are given by the dashed curves. 



A second important feature showing up in Fig. [T3] is 

that the PFA expression always overestimates the effect 
of temperature on the force between perfect (and plasma) 
mirrors. However between Drude metals, the PFA under- 
estimates this effect at small distances and overestimates 
it at large distances, the overestimation being however 
smaller than for perfect mirrors. These results clearly 
indicate that there is a strong correlation between the 
effects of plane-sphere geometry, temperature and dissi- 
pation. 

The calculation of the Casimir free energy may be done 
analytically for small frequencies corresponding to large 
plane sphere separations 



■pcrf 



3hcR^ 



AXtL^ 
cosh V A- V sinh v - 



_ 27ri 
At 

cosh V sinh^ 



(39) 



2 sinh'' V 



This simple expression is a good approximation, as 
proven by the fact that the full expression of t^p^''^ tends 
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indeed asymptotically to this simple form for small radii 
R L (dotted line on upper graph of Fig. [T21). One can 
also derive from this expression interesting information 
about the behaviour of the Casimir entropy 



-fPClf 



dT 



(40) 



This expression takes on negative values for v < 1.5, that 
is L < 1.8/im at T = 300K, which is in agreement with 
the behavior observed in the upper graph of Fig. 1131 in 
most cases d^'^'^^ decreases below unity as the distance 
increases, reaches a minimum and then increases again 
at long distances. As long as R is not too large, the 
thermal photons provide a repulsive contribution over a 
distance range that gets wider as R decreases, to become 
L < Ar/2 for very small spheres. 

We finally will compare the predictions of the 
dissipation-less plasma model and the dissipative Drude 
model for the thermal Casimir interaction in the plane- 
sphere geometry. The difference will become particularly 
clear in the high temperature limit £ 3> At where one 
only needs to take the first Matsubara frequency Co = 
when computing the Casimir free energy. In the low fre- 
quency limit, the Fresnel coefficients ([T5| for the plates 
are given by tte ~ — ''tm ~ —1 for the plasma model. 
The Mie coefficients are easily evaluated [tO, [Z1| and the 
following approximation for the Casimir force within the 
plasma model 



mcR-^ 



8XtC^ 
C ^ Xt, R, Xp , 



1 coth a 
1 + — 

a 

2nR 



Xf 



This result reproduces, as a particular case, the perfectly- 
reflecting limit when Xp ^ R. 

For the Drude model, the TE Fresnel reflection coeffi- 
cient has the well-known low-frequency limit tte ^ 0, 
whereas the TM coefficient behaves as in the plasma 
model: ttm ~ 1- The low-frequency expansion of the 
Mie coefficients are also quite different from the plasma 
case and can be found in [t^, [7i|. The resulting high- 
temperature large-distance limit for the free energy is 



Drud 



hcR^ 
'iXrC^ 



C:>Xt,R. 



(41) 



This remarkable result does not depend on the length 
scales Xp and characterizing the material response, 
whereas the corresponding plasma result (j4ip clearly de- 
pends on Xp. One can show that this is always the case 
in the high-temperature limit At ^ C. 

In the case of the Drude model with a non vanish- 
ing relaxation frequency the free energy for the Drude 
model turns out to be 2/3 of the expression for perfect 
mirrors whereas this ratio is 1/2 in the plane-plane ge- 
ometry. The latter result is explained by the fact that 
the TE reflection coefficient vanishes at zero frequency 




FIG. 14: Ratio of thermal Casimir force at T = 300/S' calcu- 
lated with the plasma model and the Drude model, as a func- 
tion of surface separation L for different radii of the sphere. 
The solid curves from bottom to top correspond to increasing 
values of sphere radii. The dashed curve is the PFA predic- 
tion. 



so that only the TM modes contribute [42, |45|. The 
change of the ratio 1/2 to 2/3 in the plane-sphere geom- 
etry has to be attributed to the redistribution of the TE 
and TM contributions into electric and magnetic spher- 
ical eigenmodes. The change is illustrated in Fig. [T4l 
where we have plotted the ratio of the thermal Casimir 
force FP^^*^ calculated with the plasma model to the one 
pTiTud obtained with the Drude model. Again, the plots 
correspond to Ap = 136nm and A^/Ap — 250. The re- 
sults of our calculations are shown by the solid curves 
with the sphere radius increasing from bottom to top. 
The ratio ^pias^^Drud yaj-jgg ^]^g plane-sphere geom- 
etry as a function of the sphere radius, clearly demon- 
strating the strong interplay between the effects of tem- 
perature, dissipation and geometry. For large spheres 
(i? ^ Ap), the ratio converges to the value 3/2, whereas 
it remains smaller for small spheres (down to 1.2 for 
R ^ lOOnm). The dashed curve gives the variation of the 
same ratio as calculated within the PFA which leads to 
a factor 2 in the limits of large distances or high temper- 
atures, corresponding to the prediction in the parallel- 
plates geometry. This factor 2 deduced within PFA is 
never approached within the calculations performed in 
the plane-sphere geometry. 



V. CONCLUSION 

In this paper we have reviewed the quantum and ther- 
mal Casimir interaction between parallel plates, corru- 
gated surfaces and plane and spherical mirrors. To per- 
form our calculations we have extensively used the scat- 
tering approach where the objects are characterized by 
scattering matrices. We have compared our results with 
predictions obtained within the PFA. When taking the 
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diffraction of the electromagnetic field correctly into ac- 
count, surprising features appear especially for the ther- 
mal Casimir force in the plane-sphere geometry, where 
the exact results differ substantially from predictions 
within the PFA. While open problems are still waiting 
to be tackled, the whole set of presented results clearly 
illustrates the usefulness and practicality of the scatter- 
ing approach in Casimir physics. 
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